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SPECIFICATION 

TITLE OF THE INVENTION 

Methods of Designing Optimal Linear Controllers 

TECHNICAL FIELD OF THE INVENTION 

This invention relates to the design of optimal multivariable linear controllers and PID 
controllers. 

BACKGROUND OF THE INVNTION 

A PID (proportional-integral-derivative) controller is a device, usually implemented in a 
computer, that is used to control a process variable (e.g., temperature, pressure, etc) of 
an industrial process. For a multiple-input and multiple-output (MIMO) process, the 
process variable y(t) and controller output u(t) are n-dimensional vector variables of the 
time t. The controller receives the sample values of y(t) and calculates its output u(t) 
according to a PID equation. The controller sends its output u(t) to the process so that the 
error e(t)=r(t)-y(t) approaches zero as time t increases, where r(t) is the set point (also 
known as the reference signal or the command signal, etc). 

There are many types of PID controllers, depending on the use of different types of PID 
equations. In discrete time form, all these PID controllers can be viewed as special cases 
of the general linear controller with the linear control equation 

Du k =Er k -Cy k 

where y k =y(t k ) is the process variable at time t k =to+kT S) to is the initial time, T s >0 is the 
constant sampling period, k=0, 1, 2, ... is a non-negative integer called discrete time 
variable, u k =u(t k ) is the controller output at time t*, and u k can be subject to upper limit 
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and/or lower limit constraints placed on one or more of its components, r k =r(t k ) is the set 
point at time t k , D, E and C are n by n matrix polynomials in the unit backward shifting 
operator z' 1 such that for any discrete time signal x kl z" 1 x k =x k .i, and one or more of the 
said D, E and C contain tuning parameters that are to be determined, 

For example, if D=(1- z" 1 )/T s *l, where I is the identity matrix of order n, and 
E=C=Ki+K p D+K d D 2 , where K», K p and K d are user specified constant coefficient matrices 
called the integral gain, the proportional gain and the derivative gain, then the linear 
controller becomes the following "type A" PID controller 



where e k =r k -y k is the error. The three gains K it K p and K d are the tuning parameters. 

If the aforesaid E is changed to Kj+K P D or Kj, then the linear controller becomes the 
following "type B" and "type C" PID controllers, respectively: 



Once its structure is selected, the performance of a linear controller or a PID controller 
depends mainly on the choice of its tuning parameters. 

How to properly choose the tuning parameters for a PID controller is a problem that has 
attracted a lot of studies ever since PID controllers became widely used in industry in the 
early 1940s. The Ziegler-Nichols tuning methods developed by Ziegler and Nichols in 
1942 and 1943 are still widely used in industry. Other model based methods choose the 
tuning parameters by minimizing some well-known control performance index such as the 



Du k = Kje k + K p De k + K<j D 2 e k 



Du k =Kie k +K p De k -K d D 2 y k 
Du k =Kie k -K p DyirK d D 2 y k 
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integral absolute errors (IAE), the integral squared errors (ISE), etc. However, practice 
shows that all these methods often lead to oscillatory control results. 

DETAILED DESCRIPTION OF THE INVENTION 

Suppose the open-loop transfer function of the process from u k to y k is A" 1 B, where A and 
B are known n by n matrix polynomials in the unit backward shifting operator z" 1 : 

A=l-A 1 z 1 -A 2 z* 2 -... -A p z p 
B= B!Z 1 + B 2 z" 2 + ... + B q z' q 

where I is the identity matrix of order n, p and q are two known positive numbers, and A^ 
A 2 , .... A p , Bi, B 2l .... and B q are constant coefficient matrices of order n. A and B can be 
obtained by use of the commercially available software product "System Identification 
Toolbox" from The MathWorks Inc. Below is another way for finding A and B. 

Select a positive integer v that is not less than n(p+q) and let k be not less than v + max{p, 
q}, where max{p, q} denotes the larger one of the p and q. Let h k be a vector defined by 

hk T =[y»<-i T yk- 2 T ... yk- P T u M T u k . 2 T ... u M T 1] 
where for any matrix X, X T denotes the transpose of X, and let 

Y=[y k y k -i... yk-vf 
H=[h k h w ... h k . v ] T 
W=diag(1, w, w 2 , w v ) 

where w is a positive real number not larger than one, and W is a diagonal matrix with 1 , 
w, w 2 , ... , and w v on its diagonal. The number k should be large enough and the controller 
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output signal u k should be chosen so as to ensure that the matrix H is of full column rank. 
Then the matrices A tl A 2> .... A p , Bi, B 2 , and B q , and the vector d can be obtained from 

[Ai A 2 ...A P B 1 B 2 ... B q d] = ((WH) + (WY)) T 

where for any matrix X, X + denotes the Moore-Penrose pseudo-inverse of X, which is a 
standard function under the name "pinv" in the software product MATLAB from The 
MathWorks Inc. (MATLAB is a registered trademark of The MathWorks Inc.). Thus the 
two polynomials A and B and a vector d are obtained. When A, B, and d are so obtained, 
the process can be approximately represented by the following linear model: 

Ay k =Bu k +d 

And the open-loop transfer function of the process is then A" 1 B. Please note that this 
invention does not deal with the problem of how to get the polynomials A and B. From 
now on it is always assumed that, one way or another, the two polynomials A and B are 
known. 

If the process is controlled by the following linear controller as described in the previous 
section: 

Du k =Er k -Cy k 

where D, E, and C are n by n polynomials in the backward shifting operator z \ and one 
or more of D, E and C contain tuning parameters that are to be specified, then the 
closed-loop transfer function from the set point r k to the process variable y k is 



Q=(A+BD' 1 C)- 1 BD" 1 E 
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From Q the characteristic polynomial in z, denoted by b(z), can be found. All roots of the 
characteristic equation b(z)=0, denoted by z<i> z 2f z 3t .... z jf where j is a positive number, 
form all the poles of Q. 

Let max {\z,l |z 2 |, |Zj|} denote the maximum of all the absolute values \z^ |z 2 |, .... and 
|zj| of the poles of Q. Then the best choice of the tuning parameters is such that the max 
{|zi|, |z 2 |, .... |Zj|} is minimized, i.e., the best tuning parameters form a solution to the 
following minimax problem: 

min max {frl, |z 2 |, .... |zj|} 

In situations where one or more of the tuning parameters must be within user specified 
regions, the above minimax problem becomes a constrained minimax problem. 

Many commercially available software products can directly be used to find the poles of a 
transfer function and the roots of a polynomial, and to solve the minimax or constrained 
minimax problem as formulated above, such as (1) the function "pole" in the "Control 
System Toolbox", which can find the poles directly from the transfer function Q, (2) the 
function "roots" in M ATA LAB, which can directly find the roots of the characteristic 
equation b(z)=0, and (3) the functions "minimax" and "fminimax" in the "Optimization 
Toolbox", which directly provides solution to the constrained or unconstrained minimax 
problem formulated above, all of which are easy to use and commercially available from 
The MathWorks Inc. 



